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Abstract 

We give sufficient conditions on initial and target measures supported 
on the sphere §" to ensure the solution to the optimal transport problem 
with the cost ^ is a diffeomorphism. 

1 Introduction 

In this paper, we will show the following two theorems. 

Theorem 1.1. Suppose that we have two probability measures fi := e^dVol and 
V := e^dVol on §" C M"+-'^, with f and g smooth functions. If 

\\Df\\ + \\Dg\\<^-^^^^ (1.1) 

TT 

where ujq satisfies 

woe"" = 2, (1.2) 

then the optimal pairing of /i and v under the optimal transport problem with 
cost given by the Euclidean distance squared on R"+^ is a diffeomorphism. 

Theorem 1.2. Suppose that we have two smooth probability measures fi := 
pdVol and v := pdVol on S" C such that 



whe- 



re 



W^|(/i, z^) < max{niin p[x), min p(x)}Ai(ri) 



Ai(n):^ . ' cos"+^0sin"-^</.d</., (1.3) 

n(n + l)(n + 2) \tt / 



'Material from all sections excluding Isection 51 represents a portion of the Ph.D. thesis of 
this author submitted at Princeton University. 
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and W2 is the Wasserstein distance, computed with respect to cost given by 
Euclidean distance squared on Then the optimal pairing of fi and v under 

the optimal transport problem is a diffeomorphism. 

Additionally, as a Corollary of Theorem 1 1.2 1 we obtain 

Corollary 1.3. Suppose that we have two smooth probability measures fi := 
pdVol andv -.^ pdVol on 5" C M"+i such that 

Hp - p||l°=(S") < niax{minp(a;), min p(a;)}A2(n) 

where 



7rVol{§"'-^)Vol{S'') f2Y+^ f"/^ „+2 



^2H^= ,\u - / cos"+^0sin"-^0d0. (1.4) 

n(n+l)(n + 2) \tt J Jq 

Then the optimal pairing 0/ /i and v under the optimal transport problem with 
cost given by the Euclidean distance squared on W^^^ is a diffeomorphism. 

For the case when the cost is given by the geodesic distance squared on S", 
regularity for general smooth positive densities was shown by Loeper in [3] . The 
cost given by the Euclidean distance squared was first investigated by Gangbo 
and McCann in [5], where the authors show examples of measures given by 
smooth densities where the optimal pairing is not given by a map. Thus there 
is a need for some condition on the two measures. 

The idea behind both proofs is essentially the same: Follow a continuity 
method, considering solutions to the elliptic optimal transport equation (de- 
fined in [section 21 below). If we can show that the image T+(a;) (also defined 
in [section 21 below) of a point a; S §" remains close enough to x, then deriva- 
tive estimates follow from arguments of Ma, Trudinger, and Wang in [5]. In 
Theorem 11.11 we obtain this closeness by first showing a gradient estimate of 
the solution u to the optimal transportation equation, while in Theorem 1 1.2 1 we 
derive the closeness directly from a bound on the T4^2~Wasserstein distance asso- 
ciated to the Euclidean distance squared on M"+^. Additionally, we are able to 
use a simple estimate to obtain a bound on the W2-Wasserstein distance from a 
bound on the L°° difference of p and p, which allows us to immediately deduce 
Corollary 11.31 from Theorem 11.21 

Finally, we present a short outline of the remainder of the paper. In lsection 21 
we give the setup of the problem. Since our cost does not satisfy the twist 
condition, we take one of the branches of the cost exponential function, then 
use that to define the elliptic equation (|2.3p . In [section 3l we calculate the 
Ma- Trudinger- Wang tensor in a specific coordinate system. In [section 4l we 
prove a gradient estimate for solutions of (|2.3p under appropriate conditions on 
/ and g. In [section 51 we show an estimate under the conditions on W2{p,p) 
given in the second theorem, and also prove a short lemma which allows us to 
obtain the estimate under conditions on \\p — p\\l°°(^s")- In Iscction 61 we show 
that the contact set for such solutions can only consist of one point, under the 
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appropriate gradient bound. Finally, in [section 71 we use the continuity method 
to prove our main theorem. 

Acknowledgement: The authors would like to thank Alessio Figalli for 
a conversation which lead to a better argument in [section 6l 

Notation: We provide here a reference table for the notation used in this 
paper. 

Notation Definition Location 



ojq Constant satisfying woe"" = 2 (|1.2p 
Ai(n) C^r' cos"+^ 0sin"- ^ m 

A.(") ^":g:;;^r"^ d )""^ rj^ cos-- ^sin-- m 

\-\ Euclidean norm on M"+^ 

(■, ■) Euclidean inner product on R"+^ 

^(•, •) Canonical metric on §" 

|-|sn Length of vectors and covectors on S" 

||Dm|| sup^gS"l-Du(x)|s.> 

.fc; (x, y) etc. ^ g|- ^c(a;, y) etc. 

d"'^ inverse matrix of c^ j 

Wtj{x) Uij{x,t) + Cij{x,T+{x)) 

Y^{x,p) Inverse of a branch of the map y i-> Dxc{x, y) = —p (|2.6p 
T+{x),T+{x) Y+{x,Duix)) before (1^ 

Y~{x,p) Inverse of a branch of the map y i— >■ Dj.c{x, y) = —p (j6.2p 



2 Set up of problem 

2.1 Monge and Kantorovich problems 

For probability measure spaces (A", ^) and (y, i^), let n(/^, be the set of prob- 
ability measures 7 on A" x 3^ such that 

^{E X 3;) = 

7(A' X i?) = 

for all C A" and E C y measurable. The cost is a measurable function 

c: X xy 

Definition 2.1. A probability measure 70 G n(/i, i^) is a Kantorovich solution 
to the optimal transportation problem between pL and v with cost c if 

/ c{x,y)d'yo{x,y) ^ inf / c{x,y)d'y{x,y). 
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Definition 2.2. A measurable map T : X ~¥ y is a Monge solution to the 
optimal transportation problem between /i and v with cost c if 

c{x,T{x))dfi{x) = inf / c{x, S{x))dii{x). 

X S#tJ.=u J X 

Definition 2.3. A real valued function u defined on X is c-convex if for each 
Xo € X , there exists a Ao G M and i/q E y such that 

u{xo) = -c(a;o,?;o) + Ao 
u{x) > — c(x, yo) + Xo, \fx E X. 

If the second inequality is strict for x ^ xo, we say the function is strictly c- 
convex. We call such a function —c{-,yo) + Aq a c-support function to u at 

Xo- 

Definition 2.4. For a c-convex function u, we define its c-transform u'^ by 

u^iy) := sup {-c{x, y) - u{x)). 
xex 

Remark 2.5. If u is c-convex, then at any fixed xo where u and c are differen- 
tiable, we can see that for some yo € y, 

Ui{xo) = -Ci{xQ,ya) 

Uijixo) > -Cij{xo,yo) (2.1) 
where the second inequality is in the sense of matrices. 

2.1.1 Twisted case 

Now suppose that for each a;o the map Dxc{xo, ■) : y ^ T*^X is defined and 
injective. In this case we can imphcitly define T{x) from a c-convex function u 
as 

-Ui{x) = Ci{x,T{x)). 
Differentiating this and taking the determinant, one obtains 

Aei{uij{x) + Cij{x,T{x))) = |det Cij(a;, r(a;))| det yQ^S.^) 
and from this we can write down the elliptic optimal transport equation: 

dei[u,,{x)+c,,{x,T(x))) = |detc,,,(x,r(a;))|-^^, (2.2) 

here /^(x) and v{y) are densities with respect to dx^ ■ ■ ■ dx" and dy^ ■ ■ ■ dy" for 
chosen coordinate systems (assuming the measures are absolutely continuous 
with respect to a coordinate volume form) . This equation is (degenerate) elliptic 
for a c-convex solution u (see Remark 12. 5p . and is invariant under a change of 
coordinates. It is a standard result that c-convex solutions of the above equation 
determine uniquely the solution to the optimal transportation problem. 
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2.2 S"- with Euclidean cost 

For the remainder of the paper, we specialize to the case 



c{x,y) := 




where |-| is the Euchdean distance on R"+^. 

In the twisted case, it is clear that each c-convex function determines a 
single- valued map T{x). This is no longer the case for this cost function on the 
sphere: For a fixed xo G S", the map Drcc{xo, •) : §" — ?• is not injective. 

In fact, if p G T*^S", \p\s^ < 1, there are exactly two points y = yi and y = y2 
such that Dxc{xo, y) = —p. With this in mind, we make the following definition. 

Definition 2.6. We define the map Y^{x,p) implicitly by 

D^c{x,Y~^{x,p)) = -p 
{x,Y+ix,p))>0 

forpeT*§", bis. < 1. 

Now for u c-convex with < 1, we use in place of the implicit 

definition of T{x) in (|2.2I) to define an elliptic equation 

detw,,{x) = |detc,,,(x,r+(x))|-^^f^, (2.3) 

where 

Wij{x) :^ Utj{x) + Cij{x,T^{x)) (2.4) 

and 

T+{x) Y+{x,Du{x)) 

which is well-defined as long as ||-Du|| < 1. Eventually, we will show that 
under appropriate conditions on the densities, solutions of this equation indeed 
determine single- valued solutions to the optimal transport problem. 

Remark 2.7. Ifu is c-convex and differentiable at Xq with \Du{xq)\§7i < 1, and 
(xq, yo) is a pair of points satisfying (|2.1[) . it is easy to see that either {xq, j/q) > 
or there exists another point yi also satisfying (j2.1[) with {xQ,yi) > 0. In the 
latter case Cij{xo,yi) > Cij{xo,yo). Hence, we see that: 

the pair {xq,T~^{xo)) satisfies (|2.ip . 
'Wij{xo) is positive semidefinite. 



3 Calculation of MTW tensor 

In this section, we utilize a coordinate system specialized to this problem on 
the sphere to calculate various quantities involving c, most notably the MTW 
tensor of [5]. 
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We will define a coordinate system centered around a point xq € S" as 
follows. First, rotate a;o so it is given by e„+i, then take coordinates for the 
upper hemisphere by representing it as the graph {x,(3{x)) over i?i(0) C M", 
where 

Note that we leave ourselves the freedom to further rotate §" as long as the 
e„+i direction remains unchanged. 

Definition 3.1. !5, Section 2] Given V,We T^E>'' and 77, C e T*E,", define 

We will say that a cost c has the property (|A3s[) at {x, y) G §" xS" with constant 
6o>0 if 

{MTW)1l ix,y)v'v^mm > '^ol^^li^l'?!!-.. (A3s) 

The following is widely known, but we carry out the calculations in our 
coordinate system here for later reference. 

Proposition 3.2. The cost satisfies (|A3sP with a uniform constant 60 = 

1 for all (x, y) e S" x §" such that {x, y) > 0. 

Proof. We will utilize the coordinate system indicated above, and calculate var- 
ious derivatives of c. First, 

_ \{x,p{x))-(.y,Hy))? 

cyx, yj — 2 

_ \x-y\^ + {Pix)~P{y)f 
2 

_ \x\^ + \y\^ - 2{x, y) + {l- \x\^ + 1 - \y? ~ 2/3 (x) /? {y)) 

2 

= l-{x,y)-l3{x)p{y). 
Thus we calculate, at generic x and y, 

Ci = ~yi - Pt {x) (3 (y) 

Ctj = -Aj (a:) /3 (y) 
Cij,k = -Pij (x) I3k (y) 

Ci,k = -Sik - Pi (x) Pk (y) 
Ci,ki = -Pi (x) Pki (y) 

Cij,ki = -Pi] (x) Pki (y) (3.1) 



and 

P^ix) 
Pij (x) 



P{x) 

'J{xj^ P{xf 
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Thus we find at a; = 0, 



Ci = -Vi 
Cij = 6ijl3 (y) 

- -A 



Ci^k ~ *^ifc 
Ci,kl = 

Ski UkVi 



Cij,kl = -Oij I ^7;^ + ^—3 . (3.2) 



Now we can calculate for V G TqS", rj G To*S", 

(MTM^)^j V'V^VkVi = -ic^,pq - c,,,rc'^^'cs,pg)cP'''c'''^VW^Vkm 



since 



> \V?\ri? 
= \V\lAri\ln 



k,l 



<i 



for any y and r]. The last equality is seen from calculation of the metric in our 
coordinates, shown below. □ 

We also make a few calculations for later use. In the above coordinates, 

9ij = Sij + Pi {x) Pj {x) 
{9ij)k = Pik (x) {x) + Pi {x) Pjk (x) 

{9ij)kl = Piki ix) Pj (x) + Pik (x) pji (x) + pii (x) Pjk {x) + Pi (x) Pjki {x) 



and at a; = we have 



9ij — Sij 

{9ij)k = 

{9ij)ki = SikSji + 6iiSjk- 
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Hence at a; = 0, 



9^' = % 
i9'')k = 

= -SikSji - Sii6jk- (3.3) 

Remark 3.3. Suppose we have a , c- convex function u. Then, when written 
in a coordinate system chosen as above centered at Xq, using the implicit relation 
in Definition \2.b\ we notice that 

4 Gradient estimate : Theorem 1.1 

In this section, we will prove an a priori gradient estimate for a solution u of our 
elliptic equation (|2.3I) . A gradient bound of the form \\Du\\ < 1 — e will ensure 
that Proposition 13.21 is applicable at (i:,T+(a:)), hence we can use the MTW 
theory to to obtain a priori second derivative estimates for u. The method we 
use is similar to that of Delanoe and Loeper in T , where the cost is the geodesic 
distance squared. In (|2.3p we will take 

v{y) :=e^'(^)^det5,,(y). 

Theorem 4.1. Let n > 2. Suppose u is a solution to equation (|2.3p such 
that 

\\Du\\ < 1. 

Then, 

Proof. Define 

where m is a solution to the equation ()2.3p . Now we let xo be the point where 4> 
achieves its maximum on and take the coordinate system defined in lscction 31 
centered at this point. We take dx^ in the direction of Du, and rotate the 
remaining n — 1 directions so that for 1 < i, j < n is diagonal at xq =. Note 
that at cc = 0, 

_ {g'"^)iUpUq + 2gP''UpUqi 

2 

= UiUii 
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and 

Vij — 2 



1 — 2 — 

p 



p 

Here we have used p.3p . Now if Ui(0) = 0, that imphes that = and u 

is constant. Thus we may assume ui(0) ^ and hence 

ui,(0) = 

for all 1 < i < n. In particular, the whole matrix Uij is diagonal at 0. 
Consider the operator 

Lv := w'-'^v^j (4.1) 

which is the second order part of the linearization of the natural logarithm 
of ()2.3p . Taking v = (j) and x = 0, and applying the maximum principle we find 
that 



> L(t){0) 

= w"-^ {-ul6iiSji +uiuiij + ^UpiUpj) 

p 

= —ulw^^ + UlW^-'uiij + uf -' UpiUpj 

p 

= -^^ + E^""("i"i-+"-)- (4.2) 

Here we have used ()3.2p and the fact that Uij is diagonal at 0. 
By differentiating the implicit relation 

Ui{x) + Ci(x,T^{x)) = 

we find that 

Thus from and Remark IH^Sl 

(T+)j(0)=u.,(0) + <S,,/3(T+(0)) 

^U^JiO)+S^J^3{Du{0)). 

In particular, (T+)^ (0) is diagonal. Note additionally, from this we find that 
"^'^ " u.j(0)+<5y/?(i?u(0))- ^^"^^ 
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Also, from p.ip we calculate 

Cyi(0,2/) = -Aji (0)/3(y) = 0. 
By the calculations at the end of lsection~3l we have 



e/(^ydetg,j(a;) 
g9(i^+(-))^detg,y(T+(a;)) 



/3(x) 



/3(T+(a;))- 



We now differentiate equation (j2.3p (after taking logarithms again) in the cci 
direction to obtain, at a: = 0, 

= w'^{u,ji + Cyi + c,j- fc(T+)J) - c'''^'(c,ij + c,jfe(T+)5^) 

- D,, if - log /3) + D,, {g o T+ - log (3 {T+) ) 

a 

+ ^ + )(T )i- 

a 

-/i + ^ + (.9ior )(r )i- 

Ml 



Substituting into (|4.2p . we obtain 



> + E + "1 ( E ^"""1 + /i - 9i{Du)l3 (Du) 



ul 



/? (i:>u) 

> V ii;"" + ui/i - uigi{Du)P {Du, 

^ P [Du) 



a>l 



Here we have used that w""(0) > and w"(0) = p(Di(oy) by ^M- Thus we 
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find that 



Hence 



{\\Df\\ + \\Dg\\)>\h\ + \g,{Du)P{Du)\ 

a>l 

> ui{n-l)iY[ w"")^ 

= (n — l)ui 

> (n — l)ui 



a>l 
^g(Du) 

min ef 



max e 



, , mine'' , ,, „ , 
(n-1) j] \\Du\\ 



\Du\\<\^]('^^] {\\Df\\ + \\Dg\\). 
n — 1 J \ mm 



□ 



5 Gradient estimate : Theorem 1.2 

We will prove the following interior estimate on Euclidean space itself and then 
show that it adapts easily to hold on embedded spheres as well. 

Theorem 5.1. Suppose that \T{x{)) — .to| = a , where T is an optimal transport 
map from p{x)dx to p{y)dy in R", with p, p supported in domains r2,r2 and 
d{xo,dft) > a. Then 

(\ l/(n+2) 
1 n(n+l)(n + 2) 1 \ 
min,gop(x) Vol{S"~^) /^/^ ^os"+2 ^sin"-^ ^ / 

where (p, p) is the Wasserstein distance between p and p. 

Proof. Rotate coordinates on E" so that xq is at the origin and T{xo) is along 
the ei axis. Define the set 

Ka = {x eW \ \x\ < acoscj)} 

where (p is the angle x makes with the ei axis. This set is a sphere of radius 
a/2 centered at the point (a/2,0, ...,0). Now take a point point xi 6 Ka- The 
monotonicity condition for optimal transport (cf. [6, Def 5.1]) says that 

(xo,T(xo)) + {xuT{xi)) > {xi,T{xo)) + {xo,T{xi)) (5.1) 
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in particular 

{xi,T{xi)) > (xi,T(xo)) = acos(j>. 
It follows that T{xi) must be in a half space and that 

\xi — T{xi)\ > acoscf) — \xi\. 

Thus we can integrate 

W2(P,P)^ / ^ p{x)dx> / p(x)dx 

> r I r^°^%(.)("-f-)%"-^d.d.siu 

Jo JS"-2 Jo 



2 

>minp(a;) ^7; I I {a cos (j) — r) r" drdasin" 

xen 2 Jq Jq 

71+2 rTT/2 

= minp(.)Fo^(S'-^)-^-^^^^ I cos"-0sin"-0.0 

and the conclusion follows. □ 

Corollary 5.2. The same estimate holds on the sphere §" C R"+^ with cost 
given by Euclidean distance squared on R"+"'^. In particular 



in II I n{n + l){n + 2) 1 2/ 

I ''ll - [mm^^npix) ^0/(8-2) eos"+2 ^sin"-^ ' '^^'^^ 



l/(n+2) 



Proof. Represent a hemisphere as a graph over the tangent space at xq. Then 
using these coordinates, repeat the above calculation. First note that the dis- 
tance of the projection in the chosen coordinates bounds from below the actual 
distance. An easy computation shows that the monotonicity condition for this 
cost (in these coordinates) is stronger than ()5.1I ) It follows that our integration 
argument over Ka is intact, noting that the volume element on the sphere is 
bounded below by dx in these coordinates. The derivative bound follows by 
noting that the norm of Du is given by the length of tangential component of 
T{x) ~x. □ 

Remark 5.3. This method can he modified to uniformly convex domains, where 
the constants depend explicitly on the upper and lower curvature bounds. The 
set Ka is less nice, but still explicit. One can repeat almost verbatim the same 
argument on a manifold with nonnegative curvature (with respect to distance 
squared cost), provided one has a lower bound on the volume element in expo- 
nential coordinates. 

We also show that a bound on the L°° distance between p and p implies a 
bound on the W2-Wasserstein distance between them, which easily allows us to 
prove Corollary [L3l from Theorem ll.il 
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Lemma 5.4. Given any two probability measures /j, — pdVol and v = pdVol 
such that p, p ^ L°°(S"), then 

W\ii,y)<TiVol (§")||p-p|U.o(s„). 

Proof. Recall by [5J Chapter 5] that 

W\ii,y)= sup -/ <f>ix)pix)dVol{x)- f i^{y)p(v)dVol{y) (5.2) 
where 

/C := {(<^, V) G C(§") X C(S") I -(/.(x) - V(2/) < c{x, y)}. 

Now note that since p and are both probability measures, we may add the 
restriction (f>{en) = tpi^n) = to the definition of the set K. without changing 
the supremum in (|5.2I) . Also, by definition we see that 



sup ( - / (t>{x)p{x)dVol{x) - I i]{y)p(y)dVol{y) 

< sup ( - / cb{x)p{x)dVol{x) - I cb^{y)p{y)dVol{y) 

0G{C(S")|0(e„)=O} 



< sup - / cp-'={x)p{x)dVol{x) - I (i}^{y)p{y)dVol{y) 

0e{C(S")|0(e„)=O} V JS" 



< sup -/ ^'^{x)p{x)dVol{x)- [~c{y,y)-cj>{y)]p{y)dVol{y) 

0G{C(S")|0(e„)=O} \ JS" JS" 

sup / (l)^'{x)p{x)dVol{x) + / cj,{y)p{y)dVol{y)] 

0G{C(S")|0(e„)=O} V JS" JS" / 

< sup (||<^^1|Lo.(s.)|b-p|U=e(s„)V0l(S")) 
0e{C(S")|0(e„)=O} 

< sup (||0^^||Lip(s„)diam(§")||p-p||io.(s„)Vol(§")) 

0e{C(S")|0(e„)=O} 

< sup p,c(-,y) II diam(§")vol(§") Hp -p|Uoe(s„) 
= 7rV^o/ (S")||p-p|Uo.(s.). 

Here we have used that since 0*^^ is a c-convex function, 
||0''||Lip(S") < sup \\D,c{-,y)\\ < 1. 

□ 



6 Nonsplitting 

First we define the c-subdifFerential of a c-convex function u at a point x by 
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Definition 6.1. 

dcu{x) := {y G §"| — c(-, y) + X is a c-support function to u at x for some A G M} . 
We also define 

T-ix):^Y-{x,Duix)) (6.1) 
where y~ is characterized by 

Dxc{x,Y^{x,p)) = -p 

{x,Y~{x,p)) <0 (6.2) 

for p e T*§", IpIs-. < 1 (compare Definition 12. 6p . Wc add the subscript u to 
emphasize the dependency on the potential function u. 

Now we show a pointwise estimate on |Z?u(a;)|s7i if dcu{x) is more than one 
point for any x. 

Lemma 6.2. Suppose that u is c-convex, C\ and \\Du\\ < ^. Then, 

dMxa) = {T+{xo)] . 

Proof. Suppose that dcu{xQ) ^ {r+(a;o)}. Then we must have T~{xo) G 
dcu{xQ). Writing y^ — T.^j^{xo), this implies that for some A G R, the func- 
tion —c{-,y^) + A is a c-support function to u at a;oj hence 

u{y~) - u{xo) > -c{y^,y^) + X- {-c{xo,y^) + A) 

= c{xo,y~) 

^ \xo -y^\^ 
2 

= 1 - cos {ds^{xo.,y^)). 

However, we also have 

u{y~) - u{xo) < ds^{xo,y^)\\Du\\ 

hence 

1 < ds-i(xo,2/")||DuI| +cos{ds'^{xo,y~))■ 

However, by the definition of T", we see that | < dsr.(a;o,y~) < tt. Thus by 
considering the real valued function ||_Du||t -I- cost on the interval tt], we see 
that 

— Iji^wll + COS — ,7r||_Du|| + cos7r| 

< 1 

by the assumption on ||-Du||, which is a contradiction. □ 
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7 Proof of Main Theorem 

By combining the appropriate gradient estimate (either Theorem 14.11 or Corol- 
lary [5?2l) with Lemma |6.2[ we can use the continuity method to show the exis- 
tence of a Monge solution to our problem, proving our two main theorems. 

Proof of Theorem \l.l[ Assume n > 2, as the case n = 1 is vacuous. 
We will apply the continuity method to the equations 

det{u,,{x) + c,,{x,T+{x))) = |detc,,,(x,T+(x))| ""^'^ ^detg^^ix) 



e9t(T+(-))^det5,,(r+(x)) 
(7.1) 



and 



det{v,,iy)+c.,iT+iy),y)) = |det c,,(T+(y), y)' e^^'^'W^^^ Mv) 



ehiTHy))^Aetg.,{T+{y)) 

(7.2) 



where, for t £ [0, 1]: 

^•'•'>-"><?^ + vw 

Let 

I :~ {t e [0, 1] I (|7.ip has a smooth solution u such that \\Du\\ < 1}. 
Now a simple calculation shows that 

WDftW < \\Df\\ 
\\Dgt\\ < \\Dg\\. 

Since e-^'dVoZgn and e^^dVols^ are probability measures, /* and gt each equal 
-log(FoZ(S")) at least once. Then, by assumption (HI]) on (||i?/|| + \\Dg\\) 
we obtain that for some < A < 1, 

min gt>~ log {Vol (§")) - A ((n - l)a;o) 
max/t < -logiVol (§")) + (1 - A) ((n - l)t^o) 
max e-^* 



< gA(n-l)a;o-|-(l-A)(n-l)cJo _ A'n.-l)uio 



min eSt 

Hence, for any t G I and a solution ut to (|7.ip . we may apply Theorem 14. II 
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to find 



\\Du,\\<[—^][^^^] {\\DM\ + \\Dg,\\) 



< 



n — 1 



TT 

_ 2 

TT 



(7.3) 



In particular, ||-Dtif|| remains uniformly bounded away from 1. Thus, from 
Proposition 13.21 and by the MTW maximum principle calculation in O Section 
4], an a priori second derivative estimate for u follows. Higher order estimates 
follow by the Evans-Krylov Theorem and standard elliptic theory, thus / is 
closed. 

To show openness, we set up the implicit function theorem as in [5J Theorem 
17.6], by taking 

G : |u e C2'"(S'") : J udVg = o| x [0, 1] |w G C"^"(S'") : J vdVg ^ o| 
to be defined as 



GKt)-e Vdctff., det(-c,,(.,T+)) ef^ ./AKW, 

At a solution u{x) at some time io: we have that G(u, io) — 0. Now the linearized 
operator on the first factor (whose principal part is a multiple of the operator L 
in (|4.ip '). is an elliptic operator with no zeroth order terms. Since the linearized 
operator has index zero, the maximum principle guarantees it is a bijection, and 
openness follows. Since = e^", we may take u = at t = and apply the 
continuity method to infer the existence of smooth solutions u to (j7.ip for all 
t e [0, 1]. Similarly, we obtain smooth solutions v to (|7.2p for all t e [0, 1]. 

We now prove the c-convexity of Ut, solutions to ()7.ip . It is clear that the 
set /' := {t G [0,1] I ut is strictly c-convex} is relatively open and contains 
0. Now take any t £ I' . By (|7.3p we may apply Lemma 16.21 to find that 
dcUt{x) = {r+(a;)} for all x £ S", that is, dcUt is a single valued map. The 
strict c-convexity of Ut implies that dcUt, hence T^^ is injective. By ()7.ip . the 
Jacobian determinant of is nonzero, so by an open-closed argument we see 
that is surjective, and hence a bijection, in fact, a diffeomorphism. Thus we 
see that for any y, 

utaT+)-Hy)) + ciiT+)-Hy),y) - 
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where {utY is the c-transform from Definition 12.41 clearly difi'erentiable by the 
above relation. Diff'erentiating this relation twice and taking determinants of 
both sides, we see that (ut)'^ satisfies equation (|7.2I) . Thus, after normalizing vt 
by adding an appropriate constant, we find that vt — (ut)*^- 

If /' ^ [0,1], let to := inf ([0, 1] \/') > 0. Then uto is c-convex but not 
strictly c-convex. By the uniform convergence of ut and vt as t — > toi we can 
see that vt„ = u^^. As above, dcU^^{y) = dcVto{y) = {T+^(y)} for all y, which 
implies in turn that ut^ is strictly c-convex, a contradiction. Thus, /' = [0, 1] 
and Ut is strictly c-convex for all t G [0, 1]. 

In particular, when t = 1, u := Ut satisfies (|2.3p . is c-convex, and dcu{x) = 
{T+(x)} for all x G Thus we may apply ^ Theorem 5.10(ii) (replacing u 
with —u)] to conclude that the measure {Id x T+)^/i is a Kantorovich solution 
to the optimal transportation problem between and v. However, since T+ is 
a diffeomorphism, we see that it is actually a Monge solution. It is unique due 
to the uniqueness of the Kantorovich solution proven in [5J Theorem 2.6]. □ 

Proof of Theorem We follow the same method as the above proof, but in- 
stead of applying Theorem 14.11 we use Corollary 15.21 and the hypotheses to 
obtain the inequality (|7.3p . The remainder of the proof follows. □ 

Proof of Corollary \1.3\ By combining Theorem ll.2l and Lemma 15.41 we immedi- 
ately obtain the claim. □ 
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